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Abstract 

We provide general methods for explicitly constructing strict Lyapunov functions for fully 
nonlinear slowly time- varying systems. Our results apply to cases where the given dynamics and 
corresponding frozen dynamics are not necessarily exponentially stable. This complements our 
previous Lyapunov function constructions for rapidly time- varying dynamics. We also explicitly 
construct input-to-state stable Lyapunov functions for slowly time- varying control systems. We 
illustrate our findings by constructing explicit Lyapunov functions for a pendulum model, an 
example from identification theory, and a perturbed friction model. 
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1 Introduction 

This paper is devoted to the study of fully nonlinear slowly time- varying systems of the form 

x = f(x,t,t/a) (1) 

for large values of the constant a > (but see Section [7] below for the extension to systems with 
controls). See Section [2] for our standing assumptions on dU). Such systems arise in a large variety 
of important engineering applications such as the control of friction and pendulums [Kj IPA21 ISlj . 
It is therefore of great interest in control engineering to develop methods for determining whether 
slowly time- varying systems are uniformly globally asymptotically stable (UGAS). When ([T|) is 
UGAS, it is also highly desirable to have general methods for constructing explicit closed form 
Lyapunov functions for See for example [S3J lASWj |M| IMM3j IMMDj for discussions on 
the essentialness of Lyapunov functions for feedback design and robustness analysis. See also 
[MMD[ IPA1] for the dual problem of stabilizing rapidly time-varying systems, and see Remark [8] 
below for the relationship between our methods for constructing Lyapunov functions for rapidly 
and slowly time- varying systems. 

One popular approach to studying (pQ) is to first establish exponential stability of the corre- 
sponding "frozen dynamics" 

x = f(x,t,r) (2) 
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for all relevant values of the parameter r including cases where the exponent in the exponential decay 
estimate can be negative or positive for some values of r but is positive on average [K ^ IPA2|,IS1| . The 
stability of the frozen dynamics is then used to establish stability of (Q}. However, these earlier 
results do not lead to explicit Lyapunov functions for ([1]) that would be needed for robustness 
analysis. The main goals of our work are (i) to show that explicit Lyapunov functions for (Q} 
can be explicitly constructed in terms of a suitable class of oftentimes readily available Lyapunov 
functions for ([2]) when the constant a > is large enough, and (ii) to show how to relax the 
exponential like stability assumptions on ([2]) and also allow cases where r is a vector, thereby 
broadening the class of dynamics to which the frozen dynamics method can be applied. 

The rest of this paper is organized as follows. In Section [2] we provide the relevant definitions 
and standing assumptions on (P) . We state and prove our main result in Sections [3HH In Section 
we extend our main result to cases where the Lyapunov functions for the frozen systems satisfy less 
restrictive properties than those in Sections OH In Section [61 we illustrate the wide applicability 
of our results using four examples. In the first two examples, the family of Lyapunov like functions 
for ([2]) are independent of r, so the strict Lyapunov functions we construct for (pQ) are valid for all 
a > i.e. (pD) is UGAS for all a > 0. Our next two examples involve a mass spring model with 
slowly time-varying coefficients from [DDNZ and an identification model similar to those studied 
in [PAlj . and illustrate the more general situation where ([T]) is not necessarily UGAS for all values 
of a > 0. In each case, the dynamics have slowly time-varying coefficients and so are beyond the 
scope of the previously known Lyapunov construction methods. In Section [TJ we show how to 
extend our results to systems with controls using input-to-state stability. We close in Section [8] 
with some suggestions for further research. 

2 Definitions, Assumptions, and Lemmas 

We let /Cqo denote the set of all continuous functions p : [0, oo) — ► [0, oo) for which (i) p(0) = and 
(ii) p is strictly increasing and unbounded. We let ICC denote the class of all continuous functions 
(3 : [0, oo) x [0, oo) —* [0, oo) for which 

(I) /3(-,t) G /Coo for each t > 0, 
(II) f3(s, •) is non-increasing for each s > 0, and 
(III) f3(s, t) -> as t -> +oo for each s > 0. 

When we say that a function p is smooth (a.k.a. C 1 ), we mean it is continuously differentiable, 
written p G C ■ (For functions p defined on [0, oo), we interpret p'(0) as a one-sided derivative, 
and continuity of p' at as one-sided continuity.) We let | • | denote the Euclidean norm. A 
continuous function p : [0, oo) — > [0, oo) is called positive definite provided it is zero only at zero. 
When r i— > p(r) G M rf is a function with differentiable components, we use p'(r) to denote the vector 
(p[(r),....,p' d (r)). 

The following definitions and lemma apply to general nonlinear systems 

x = h(x, t) (3) 

evolving on the state space W 1 where h is locally Lipschitz (but see Section [7] for the extension to 
control systems). Later we specialize to systems with multiple time scales and frozen parameters, 
e.g. h(x,t) = f(x,t,p(t/a)) or h(x,t) = f(x,t,r) for given constant parameters a and r and 
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suitable functions p. We always assume ([3]) is forward complete meaning for each x Q € R n and 
t Q £ IR>o := [0,oo) there exists a unique trajectory 

[t , 00) B t 1 ^ (j)(t;t ,x ) 

for © that satisfies x(i G ) = x D . We assume all of our uncontrolled dynamics (J3J) are uniformly 
state bounded meaning there exists o^ € /Coo such that \h(x,t)\ < ah(\x\) everywhere. 

Definition 1 We say that |3p is uniformly globally asymptotically stable (UGAS) provided there 
exists P E fCC such that 

\<P(t;t ,x )\<p(\x \,t-t ) (UGAS) 
for all x £ M n , i D € [0, 00), and t > t Q . 

Definition 2 A smooth function W : M. n x [0, 00) — > [0, 00) is called a Lyapunov function for |3|) 

provided there are functions 01,02 £ ^00 ft positive definite function 03 smc/i that 



(LI) ai{\x\) <W{x,t) <a 2 {\x\) and 
(L2) W t {x,t) + W x (x,t)h{x,t) < -o 3 (M) 



hold for allt>0 and x G R n . 



The subscripts on W denote partial gradients. In what follows, we often omit the arguments 
x, t, etc. in our functions when they are clear from the context; and all (in)equalities should 
be interpreted to hold wherever they make sense. A smooth function W : W 1 x [0, 00) — > [0, 00) 
that admits 01,02 £ K-oc such that (LI) holds everywhere is called uniformly proper and positive 
definite. The following lemma is standard |ELWl IK] : 

Lemma 3 // (0j admits a Lyapunov function, then it is UGAS. 

A simple application of Fubini's Theorem yields the formula 

/ f e(l)dlds= f (r - t + c)@(r)dr 

Jt—c J s J t—c 

and therefore also the following [MM3]: 

Lemma 4 Let B : M — > M be continuous and bounded in norm by some constant M > 0, and c > 
be given. Then 



ft i-t 



(A) 



(B) 



Q(l)dlds 



d 



t—c J s 

t ft 



c 2 M 
< and 



dt J t 
hold for all t £ 



0(l)dlds = cS(t) - / 9(r)dr 



c J s 



t-c 
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3 Statement of Main Result and Remarks 

For simplicity, we assume our system ([1]) has the form 

x = f{x,t,p{t/a)) (4) 

where p : R — » R d (for some integer d) is bounded and its components Pi,---,Pd have bounded first 
derivatives. We set 

p := sup{|p'(r)| :r£R) and K(p) := {p(t) : t £ R} . 

Our next assumption is a variant of those of [PA2J Theorem 2] (but see Section [5] for results under 
weaker assumptions). 

Assumption 5 There exist ai,a 2 G fCoo, positive constants c a , q>, and T, a continuous function 
q : R d -> R, and a C 1 function V : R n x [0, oo) x R d -> [0, oo) such that 

A! ai (\x\) < V(x,t,r) < a 2 {\x\), 

A 2 V t (x, t, t) + V x (x, t, r)f(x, t, t) < -q(r)V(x, t, t), 

^3 \V T (x,t,T)\ < c a V(x,t,T), and 

Aa ft_ T q{p(s))ds > c b 

hold for all x G R n , t > 0, and r G TZ(p). 

Note that A 2 is weaker than the standard exponential stability property of the frozen dynamics 
since we do not require at or a 2 to be quadratic functions and moreover g(r) can take non-positive 
values for some choices of the vector parameter r. However, A4 requires that q be positive on 
average along the vector p(s). 

Theorem 6 If satisfies Assumption^ then for each constant a > 2Tc a p/c\ ) , the dynamics 
are UGAS and 

f— c— 

T I I q{p{l))dlds 

V*(t,x):=e J i~ TJs V(x,t,p(t/a)) (5) 



is 



a Lyapunov function for Q). 



Remark 7 Compared to the known results JK\ \PAfy \S1]/ , the novelty of Theorem is that (a) 
we allow fully nonlinear systems including cases where the function q can take both positive and 
negative values (which corresponds to the allowance in lSl\j of eigenvalues that wander into the right 
half plane while remaining in the strict left half plane on average) and (b) we provide an explicit 
Lyapunov function ^ for the original slowly time-varying dynamics. In general, the conclusion of 
Theorem^ may or may not hold for small values of a. We illustrate this in Section^ below. 

Remark 8 A completely different approach to slowly time-varying systems x = f(x,t,t/a) (for 
large constants a > 0) is to transform the system into a rapidly time-varying system and to then 
try to construct a Lyapunov function for the resulting rapidly time-varying system directly. The 
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transformation is done by simply setting s = t/a which gives rise to the new rapidly time-varying 
system 

x(s) = g(x(s), s, as) := af(x(s), set, s) (6) 

in terms of the new reseated time variable s. However, it might be difficult to apply the Lyapunov 
function construction methods of JMMDlj or other known methods to build an explicit Lyapunov 
function for {6]j. This is because these earlier results are for fast time-varying dynamics having a 
different form from |S]) and moreover they require Lyapunov functions for so-called limiting dynam- 
ics; cf. \PA 11 Property 2], and see \MMLA Section 3.1] for a generalization in the same vein. This 
motivates our direct construction of Lyapunov functions for slowly time-varying dynamics, which 
may be viewed as a complementary approach to the time rescaling method since we do not require 
limiting dynamics. 



4 Proof of Theorem [6] 

By A2-A3 and our choice of p, the time derivative of 

V(x,t) := V(x,t,p(t/a)) 

along the trajectories of (H|) satisfies: 

V = Vt(x,t,p(t/a))+V x (x,t,p(t/a))f(x,t,p(t/a))+V T (x,t,p(t/a)) 



(7) 



p'(t/a) 



< -q(p(t/a))V(t,x) + V T (x,t,p(t/a)) 



p'(t/a) 



a 



< 



-q(p(t/a)) + V(x,t). 
a 



a 



To simplify the notation, let us define 

E(t, a) := e 
Since p is bounded and q is continuous, 



q{p{l))dl 



ds 



S(t) := q(p(t)) 



(8) 



(9) 



is bounded in norm by some constant M > 0. Therefore, along the trajectories of (j4]), Lemma [4] 
(B) with the choices ([9]) and c = T gives 



V* = E(t,a) 



V + q(p(t/a)) - - J^q(p(l))dl | V 



Substituting the formula for V, it follows from A4 that 

Vl < E(t,a 
< E(t,a 



CaP _± A 

a T a 

Cap C b 

La T 



Jl_ T q{p(l))dl 



V 



(10) 



V(x,t). 
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Applying Lemma 0] (A) with the choices Q and c = T gives 

e a ™/ 2 >£(t,a)>e- a ™/ 2 
everywhere. Hence, for a > 2Tc a p/c} 1 , (fTUj) gives 

V2(x,t) < -A e -^*^( x , t )< -asdxl), (11) 

where 03 (s) = ^e~ a ™ / 2 a\{s) is positive definite; and 

&i(M) < Vl(x,t) < a 2 (\x\) (12) 

everywhere, where 

&1 ( s ) := e~ a ™l 2 oi{s), a 2 (s) := e" TA? / 2 a 2 ( S ) 

are of class /Coo. Since (fTTj) - (fT2j) imply that Va is a Lyapunov function for fl3J), Theorem [6] follows 
from Lemma EJ 

5 More General Families of Lyapunov Functions 

We next show how to relax requirements A2-A3 from Theorem [6l We continue to use the notation 
we introduced in Section [3J We assume the following in the rest of this section: 

Assumption 9 There exist 0,1,0,2 G /Coo, a positive definite C 1 function fi, positive constants T, 
c a , and db, a continuous function q : R rf — > R, and a C l function V : R n x [0, 00) x R d — > [0, 00) 
such that 

f r 1 

lim / -—dl = +00 (13) 

and 

M oi(\x\) < V(x,t,r) < 5 2 (|x|), 

1 2 T4(x, i, r) + K(x, i, r)/(x, t, r) < -g(r)/i(F(x, t, r)), 

^L 3 |K(x,t,r)| < c a /i(V'(x,i,r)) J and 

^4 ft_ T q(p(s))ds > c b 

hold for all x G R n , i > 0, and t e Tl{p). 

Notice that Assumption is the special case of Assumption in which fj,(l) = I. We prove the 
following: 

Theorem 10 If ^ satisfies Assumption^ then there exists k £ /Coo for which the requirements 
of Assumption^ are satisfied with V := k(V). Therefore, for each sufficiently large value of the 
constant a > 0, the dynamics ^ is UGAS and admits a Lyapunov function of the form 

It suffices to prove the first statement of Theorem [10] since the second statement is immediate 
from Theorem [6) To this end, we use the following important observation: 
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Lemma 11 If fi £ C 1 is positive definite, then 

lim I -^—dl = -oo. (14) 

Proof: Since /i£C' and /x(0) = 0, we can find C3 > such that fj,(r) < c^r for all r G [0, 1]. Hence, 
for each r £ (0, 1], we get 

1 1 f 1 1 1 

dl > —dl = ln(r) (15) 



It follows that, for all r £ (0, 1], 



S'wf £ ^< r >' < 16 > 



Since lim r _ > o+ ln(r) = —00, the lemma follows. 

Given a constant £ > which we specify later, it follows from (|13p and (|14p that 



1 -dl 



e A M(0 , r > (17) 
0, r = 



is continuous and unbounded. In particular, k(r) — * as r — > + . Set B = sup{^'(s) : < s < 1}, 
which is positive since [i is positive definite. 



Lemma 12 The function |j?[ ) is C 1 u>/ien £ = 25. 
Proof: It suffices to prove that 

k'(r) as < r 0+, (18) 
since k(0) = 0, because then k(r)/r — > = A;'(0) as r — ► + . To this end, first note that 

C 1 d< 

fc'(r) = -4-re Jr K l ) , Vr > (19) 
H[r) 

and that for all r G (0,1], we have 

5HMD)-N/.W)] = / 1 ||<«</ 1 ^ 



by our choice of B. Since /i is positive definite and £ is positive, this implies 

r 1 di 

_J_ e ^Pn0*(l))-ln(/*(r))] > J^ e '^ J r W) = k! \r) 

H(r) n(r) 

for all r £ (0, 1], i.e., 

^(l))^//^)* -1 > fc'( r ) > Vr 6(0,1] (20) 
so (|18p follows from our choice of £. ■ 
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The fact that k 6 /Coo n C 1 is now immediate from Assumption [9] and Lemmas [IT] and [I2j Let 
us now verify that k satisfies the requirements of Theorem [10J From the definition of k, we deduce 
that 

r 1 i 

9R -2B / — ——dl 
k'(V)fi(V) = -=£-e Jv W) fi(y) = 2Bk(V) 

n(y) 

when V 7^ 0. Therefore, by assumption A2, the time derivative of V = k(V) along the trajectories 
of ([!]) satisfies 

V t (x, t, t) + V x (x, t, r)f(x, t, t) = k'{V{x, t, r))[V t (x, t, r) + V x (x, t, r)f(x, t, r)] 

< -k'(V(x,t,T))q(T)fl(V(x,t,T)) 

= -2Bq(r)V(x,t,T) 

everywhere, and condition A3 from Assumption implies 

|K| = k'{V)\V T \ < c a k'(V)fi(V) = 2Bc a V 
everywhere. Therefore Assumption [5] holds using 

Oi(s) := koai(s) for % = 1,2, V := k(V), q(r) := 2Bq(r), c a :=2Bc a , c b :=2Bc b . (21) 
The result now follows from Theorem [6l 

6 Examples 

We illustrate our constructions using four examples. In the first two examples, the functions V 
from Assumption [5] do not depend on the frozen parameter r, so we can conclude that §S§ is UGAS 
for all values of the constant a > 0. We then turn to a slowly time- varying friction dynamics 
and an example from identification where V depends on r, and where we can consequently only 
conclude the UGAS property of (jU) when a > is sufficiently large. Set V(x,t,r) := Vt(x,t,r) + 
V x (x,t,T)f(x,t,T) everywhere. 

6.1 Stability for all a > 0: a scalar example 

Consider the one-dimensional system 



y/1 + x 2 



I - 90 cos 2 ( - 
a 



(22) 



and the uniformly proper and positive definite function 



V{x,t,T) = V(x):=e VT +^-e. (23) 
Let us verify Assumption [5] for this choice of V and the frozen dynamics 

x = f(x,t,r) := _£_[1-90t], 0<r<l. (24) 
yl + r 
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This gives 



V(x,t,r) 



1 + x 2 



[1 - 90r] 



„2 
1 + X 2 



(25) 



Simple calculus calculations everywhere give 



2e^ 



jV(x) > Y^e^ > -V{x) 
2e^ 



so (|25p everywhere gives 



V(x,i,r) < 



e- 1 



45r 



y(x). 



(26) 



(27) 



Moreover, for each t > 0, we get 



t— 7T 



45 cos 2 (s) 



2e^ 
e- 1 



/45 2e^\ 

T " ^1 > 



which shows that Assumption [5] is satisfied. We conclude from Theorem [6] that for large enough 
constants a > 0, (|22p is UGAS and has the Lyapunov function 



45 cos 2 (Z) 



2e 



>/2 



e- 1 



ds 
J V"(i) 



-^(f)-Hr-S 



45(e-l) 



(28) 



where V is in (|23p . In fact, since in this case V does not depend on r, it follows from our proof 
of Theorem [6] that for any constant a > 0, the system (f22j) is UGAS and admits the Lyapunov 
function (l28l). 



Remark 13 T/ie dynamics A22\) illustrates the fact that our approach applies to systems which are 
not globally exponentially stable. Indeed, it is clear that (2<fy is not globally exponentially stable 
since its vector field is bounded in norm by the constant 91. 

6.2 Stability for all a > 0: a pendulum example 

Our constructions also apply to the slowly time- varying pendulum dynamics [PA2] 



Xi 

i'2 



X2 



-xi - [1 + b 2 (t I a)m(x , t)]x2 



(29) 



assuming 
(VI) m : R 2 x 



[0, 1] is Lipschitz continuous; and 



(V2) 62 : K — ^ (—oo,0] is globally bounded, and there are positive constants T and q, such that 
5 + T j£_ T 6 2 (0 > c b for all t G E. 

The dynamics ()29p was shown to be UGAS for certain choices of the function 62 in [PA1] ; see [K] 
for related results, and |Slj for results that are restricted to the linear case. However, these earlier 
results do not lead to explicit Lyapunov functions for ([29|) . In order to build Lyapunov functions 
for (|29p for large constants a > 0, we use the following observation for the corresponding frozen 
dynamics f(x,t,r) := (x2,—x\ — [1 + rm(i, t)]x2)' 
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Lemma 14 The function V (x) := x\ + x\ + x\X2 satisfies W(x)f(x,t,r) < — [1 + 5t]V(x) for all 
x £R n ,t£R, and r < 0. 

Proof: By grouping terms, one easily shows that 

W(x)f(x, t, t) = — V(x) — 2rm(x, i)^ — rm(x, t)x\X2 (30) 

everywhere. Since 

V > x\ + x\ - \x\X2\ > -^x\ + -x\ > \xix 2 \ 

everywhere, we get 

— 2Tm(x,t)x 2 < —4rm(x,t)V(x), —Tm(x,t)xiX2 < —Tm(x,t)V(x) (31) 

everywhere. The lemma follows by substituting (|3ip into (|30p and recalling that < m(x, t) < 1 
everywhere. ■ 
The following is an immediate consequence of Lemma [T4"l the proof of Theorem [6l and the fact 
that V(x) := x\ + x\ + x\X2 only depends on x: 

Theorem 15 Let f23\) satisfy (Vl)-(V2). Then $2§\) has the Lyapunov function 

t_ t_ 
^ j f a b 2 {l)dlds 
Vl{t,x):=e J i~ TJs {x\ + x\ + Xl x 2 ) (32) 

for each choice of the constant a > 0. Hence, \29\) is V GAS for all constants a > 0. 
6.3 Friction example revisited 

We next illustrate Theorem [6] using the one degree-of- freedom mass-spring system jDDNZl IMMD] . 
The mass-spring system arises in the control of mechanical systems in the presence of friction. 
However, in contrast to [MMD] where the dynamics are assumed to be rapidly time- varying, here we 
consider the case where the dynamics are slowly time- varying. While slowly time-varying dynamics 
can be transformed into rapidly time-varying dynamics by rescaling time, doing so for the slowly 
time-varying mass spring system produces a new dynamic that does not lend itself to the known 
methods; see Remark [8] above for details. For this reason, we directly apply the slowly time- varying 
theory we developed in the preceding sections. 

Let us recall the model [MMD] . The dynamics are given by 

XI = X 2 . , 

x 2 = -a 1 (t/a)x 2 -k(t)x 1 -{a2(t/a)+a 3 (t/a)e~- (3l ^}sa,t(x2) V ; 

where x\ and X2 are the mass position and velocity, respectively; <7j : [0, oo) — > (0,1], i = 1,2,3 
denote positive time-varying viscous, Coulomb, and static friction-related coefficients, respectively; 
j3\ is a positive constant corresponding to the Stribeck effect; /i(-) is a positive definite function also 
related to the Stribeck effect; k denotes a positive time-varying spring stiffness-related coefficient; 
and sat(-) denotes any continuous function having these properties: 

(a) sat(0) = 0, (b) £sat(£) > G R, 

(c) lim sat(£) = +1, (d) lim sat(£) = -1 ^ 

£— »+oo £— > — oo 
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Following [MMDJ, we model the saturation differentiably as 



sat (#2) = tanh(/?2X 2 ), (35) 

where j3 2 is a large positive constant, so |sat(x 2 )| < /? 2 |x 2 | for all X2 £ K. However, unlike |MMD| . we 
assume the friction coefficients vary in time slower than the spring stiffness coefficient so we restrict 
to cases where a > 1. We are going to establish the stability of (I33p and construct corresponding 
Lyapunov functions V a when the constant a > is sufficiently large. 

Our precise mathematical assumptions on ([33]) are: k and the cjj's are C 1 functions with bounded 
derivatives; /i has a globally bounded derivative; and there exist constants Q,,T > such that 

t 

oi{r)dr > c b Vt > 0. (36) 

t-T 

We also assume this additional condition whose physical interpretation is that the spring stiffness 
is nonincreasing: 

3k Q ,k > s.t. k a < k(t) < k and k'(t) < Vt > 0. 
The frozen dynamics x = f(x,t,T) for (J33J) are 

xi = x 2 

x 2 = -T lX 2 - k(t) Xl - {t 2 + T3 e~M x ^} sat(x 2 ) 

where r = (ti,T2, T3) 6 [0, l] 3 is now a vector of parameters. We apply our construction from 
Theorem [6] with p(t) = <r 2 (t), 03(£)) and the function 

F(x,i,r) = A(yfc(i)xf + x|) + rix 1 x 2 where A= 1 + ^+ (38) 

We first verify the conditions of Assumption O Since A > max{l, l/k } and n < 1, we have 

^(x 2 + x?,) < F(x,t,r) < A 2 k(\xi\ + |x 2 |) 2 < 2yl 2 A;|x| 2 (39) 

everywhere. Let us now compute V(x,t,r) for all values r E [0, l] 3 . Since fc'(t) < everywhere, 
this gives 

F(x,i,r) < K.(x,i,r)/(x,i,r) 

= [2Ak(t)x! + rix 2 ]x 2 - [2Ax 2 + rixi]{nx 2 + [r 2 + T ^-^ x ^] sat(x 2 ) + fc(t)xi}. 

Therefore, by grouping and canceling terms, we also have 



V(x,t,r) < -nkoxj - (2At! - Ti)xl + n(l + 2f3 2 )\x lX2 \ 

n%x 2 + (A - i/2)nx 2 - n(l + 20 2 )\x 1 x 2 \ 



h 



feo, ,2 fco /, , l + 2/?2, V, /n(l + 2/? 2 ) 2 TJ \ a 

= " T1 T N " ri T V |xi| "^^ k2| J H — 2£ — + T -M)* 2 

< -Bl^^ 

where the first inequality follows from (|34|) (b) . the inequality | sat (2:2) | < /5 2 |x 2 |, and the fact that 
Tj 6 [0, 1] for each i; the second inequality uses the fact that A — 1/2 > & Q /2; and the last inequality 
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is from the choice of A and the bounds (|39p. Hence, Assumption [5] of Theorem [6] readily follows 
from (|36p with the choices 

q{T) = AJVk> Ca = l - 

We conclude as follows: 

Corollary 16 Under the preceding assumptions, there exists a constant a > such that for all 
constants a > a a , the system 133\) is UGAS and admits the Lyapunov function 



T 



-T Js 



a±(l)dl ds (40) 



V a (t, x) := V(x, t,p(t/a)) e 
where V is the function defined in fc38\) . b = k /(AA 2 k), and p(t) = (o~i(t), 02(f), a^t)). 



6.4 Identification dynamics revisited 

Our Lyapunov function constructions also apply to the slowly time-varying dynamics 

x = h(t/a)m(t)m T (t)x, i£l n (41) 
assuming there are positive constants T, c, a, and a such that 

(XI) h : R — ► [—a, 0] is continuous with a bounded first derivative and f._ T h(r)dr < —a for all 
t e R. 

(12) m : R — > M n is continuous and satisfies |m(t)| = 1 and al < j/ +C m(r)m? (r)dr < al for all 
i e R. 

where I is the identity matrix and for matrices A,B E M nxn we use B — A > to mean that B — A 
is positive semi-definite. See also Remark [23] below for the generalization of our result to control 
affine systems 

x = /i(i/a)m(f)m T (t)a; + g(x,t,t/a)u 

for suitable matrix valued functions g. The particular case x = —m(t)m T (t)x of (JHJ) has been 
studied extensively in the context of identification theory |MMD[ PAl] . In [MMD] , we showed how 
to construct explicit Lyapunov functions for the rapidly time- varying system x = f (at)m(t)m T (t)x 
for appropriate nonpositive functions / and large positive constants a. However, these earlier 
results do not lead to explicit Lyapunov functions for the slowly time-varying dynamics (|41j) for 
large constants a > 0; see Remark [8] above. Instead, we construct Lyapunov functions for ([IT]) 
using the following analogue of [MMD] Lemma 6]: 

Lemma 17 Assume there are positive constants T , c, a, and a such that (II) -(12) are satisfied 
and set 

ft rt g Q,2g4 

P(t,r) = Kl-T / m(l)m T (l)dlds, where k = — I he 2 . (42) 

Jt-c Js 2 2 « 

Then for each t £ [—a, 0] the function 

V(x,t,r) = x T P(t,r)x 

satisfies the requirements of Assumption^ for the frozen dynamics f(x,t,r) = Tm(t)m T (t)x and 
p(s) = h(s). 
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Proof: We first apply LemmaU^B) to 0(i) = x T m(i)m T (i)a; for each x and r G [— a, 0] to get 



(x, t, r) = — tcx t m(t)m T (t)x + rx 



m(l)m J (l)dl 



t-c 



X 



and 



<9x 



(x,t,T)f(x,t,T) = 2TX 1 



/, / - r / / m(l)m T '(/) dZ ds 
t— 5 •/ s 



m(t)mJ (t)x 



so the derivative V" = ^ + ^/ along the trajectories of / satisfies 

V" < t[(2k — c)\m T (t)x\ 2 + a|x| 2 ] + T 2 \x\\m T (t)x\c 2 
= t{(2k — c)\m T (t)x\ 2 + a\x\ 2 + T\x\\m T (t)x\c 2 } 



(43) 



where the inequality follows from LemmaS^A), (12), and the facts that |m(t)| = 1 and r < 0. Set 

a 



2c 2 



a 



Then the triangle inequality gives 



1 



i T/\iii i 1 9 i T~ / \ 1 1 

to (t)x \\x\ < uj\x\ H m (t)x\ 

4u) 



so since r < 0, we get 



(2k — c)\m 1 (t)x\ 2 + a\x\ 2 + t\x\ \ m T (t)x\c 2 



> (2k - c)\m T (t)x\ 2 + a\x\ 2 + ujtc 2 \x\ 2 + ^—\m T (t)x\ 2 



> 



2k- c + 



a\x\ 



\m T (t)x\ 2 + (a + wrc 2 )|x| 



by our choices of k and u. This and (|43l) give V < -f=|x| everywhere; and LemmaUJA), the fact 
that = 1, and our choice of k give 



k\x\ 2 < V(x,t,r) < (k + c 2 o)\x\ 

dV, , - 1 
Q^(x,t,r) 



< ~\X\ 2 <V(x,t,T). 



Hence, 



V(x,t,r) < —q(r)V(x,t,T), where q(r) :- 



to 



2(k + c 2 a) 

everywhere. Therefore, we can satisfy the requirements of Assumption [5] using p(s) = h(s) and 

„2 



ai(s) = ks , ct2(s) = (k + c a)s , c a = 1, c& 



a 



2(k + c 2 a) 



(44) 



This proves the lemma. 

The following is an immediate consequence of the preceding lemma and Theorem [6) 
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Theorem 18 Let \J^1\ ) admit positive constants T, c, a, and a such that (Z1)-(X2) are satisfied 
and choose c& as in Then for any constant a > 2Tsup{|/i'(r)| : r G R}/q,, the function 



2T(re + c^(i) 



-T Js 



h(l)dl ds 



V(x,t,h{t/a)) 



(45) 



with V(x,t,r) = x T P(t,r)x as defined in Lemma 17 is a Lyapunov function for fijjty - Hence, (L/i[ ) 
is U GAS for all constants a > 2Tsup{|/i'(r)| : r G R}/c&. 



7 Input-to- State Stability 

We next extend our results to control affine slowly time-varying systems 

x = f{x,t,p(t/a)) + g(x,t,p(t/a))u (46) 

evolving on R n with control values u G R m , where / : M. n x [0, oo) x M. d — » R ra and g : M. n x [0, oo) x 
M. d — > R nxm are locally Lipschitz functions that admit 04 G /Cqo such that 

|/(ar,i,p(t/a))| + \g(x,t,p(t/a))\ < a 4 (|z|) 

everywhere, and where p : R — > M d for some d is bounded with a bounded first derivative. The 
control functions (i.e. inputs) for (|46[) comprise the set U of all measurable essentially bounded 
functions u : [0, 00) — ► R m with the essential supremum norm | • |oo- We assume throughout 
this section that Assumption [5] holds for some function V G C 1 and that ai G /Coo and c a from 
Assumption [5] are also such that 

A 5 \V x (x,t,p(t/a))\ < c a y/ai(\x\) 

A 6 \g(x,t,p(t/a))\ < c a (l + </ai(|x|)) 

hold for all t > 0, a > 0, and a; G R n . Notice that A5 reduces to a linear growth condition when 
«i(x) = \x\ 2 and so automatically holds in the classical case where V has the form x T P(t)x for a 
suitable bounded positive definite matrix. We show that when Assumption [5] and A§-A§ hold, and 
when the constant a > is sufficiently large, the control system (|46p satisfies the input-to-state 
stable (ISS) property and admits the ISS Lyapunov function ([5]). We first recall the relevant ISS 
definitions from [ELWl lMM2l El [S3]. 

For a general locally Lipschitz control affine system 

x = h(x,t) + J(x,t)u (47) 

where h and J are uniformly state bounded (as defined in Section [2]), and for given values t > 0, 
x G R n , and u G IA, we let i 1— > 0(t;t o ,x o ,u) denote the unique maximal solution of the initial 
value problem 

x(i) = h(x(t), t) + J(x(t), t)u(t) a.e. t, x(t ) = x$ . 

We always assume all trajectories (p(-;t ,x ,u) so defined are defined on all of [t a ,oo). Later we 
specialize to the controlled system (146 p for fixed constants a > 0. 

Definition 19 We say that f^?] ) is input-to-state stable (ISS) provided there exist (3 G ICC and 
7 G /Cqo such that 

\<f>{t;t 0t x ,u)\ < /3(|x |,i-t ) + 7 (|u|oo) (ISS) 
holds for all t > t Q , t Q > 0, x Q G 1", and u G IA. 
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Definition 20 A smooth function W : W 1 x [0, do) — > [0, oo) is called an ISS Lyapunov function 
for |7?p provided there exist functions a%, a^, X £ ^-oo and a positive definite function Q3 such that 

(1) ai{\x\) < W(x,t) < a 2 (\x\) and 

(2) \u\<x{\x\) implies W t (x,t) + W x (x,t)[h(x,t) + J(x,t)u] < -a 3 (\x\). 
hold for all x G W n , t>0, andueR" 1 . 

The following lemma comes from [ELW| : 

Lemma 21 If J^7| j has an ISS Lyapunov function, then it is ISS. 

We prove the following analogue of [MMDj Theorem 4]: 

Theorem 22 Assume satisfies Assumption\M Assume A^-Aq everywhere hold where c a , V, 
and a\ are chosen as in Assumption^ Then for each constant a > 4Tc a p/c b , the dynamics 
are ISS and 



V*(t,x) :=e 
is an ISS Lyapunov function for UM. 



-T Js 



q(p(l))dl ds 



V(x,t,p(t/a)) 



(48) 



Proof: We indicate the changes needed in the proof of Theorem [H Consider the function x £ 
defined by 



c h ^ai(s) 



2Tcl(l+^{s) 



X(s) 



where a\ and c a are as in Assumption El This function is of class /Coo since s 1— > a^ 4 (s) and 
s 1 ^ s 2 /(l + s) are both Koo- Our assumptions imply that if \u\ < x(M)> then 



\V x (x,t,p(t/a))g(x,t,p(t/a)u\ < 



c b ai(\x\) 



2T 



Cb 



< ^V(x,t,p(t/a)) 

everywhere. Define V as in ([7]) and E(t,a) as in (JS|) . Then, along any trajectory x = (j){t) of (j4U 
with inputs u satisfying | 

u|oo — x(|0(^)I) everywhere, wg get 



v 



j(p(t/a))V(t,x) + 



C a p C b 

a 2T 



V(x,t) 



everywhere and therefore (by reasoning exactly as in Section S]) also 



Vl < E(t,a) 
< E(t,a) 



c a p c b 1 
a 2T T 



t/a 



t/a-T 



j(p(l))dl 



V(x,t) 



CgP C b 

a 2T 



y c b E(t, a) y 
V(x,t) < — — V(x,t) 



when a > 4Tc a p/c b , where the second inequality follows from A4, and the last inequality follows 
from our choice of a. We then argue exactly as before to show that is an ISS Lyapunov function 
for (|46p for all constants a > ATc a p/c b . The theorem now follows from Lemma [2T1 I 
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Remark 23 Theorem [2H readily applies to our friction example from Section \6.3\ above, show- 
ing that is actually an ISS Lyapunov function for the slowly time-varying controlled friction 
dynamic 

x 2 = -a 1 (t/a)x 2 -k(t)x 1 +g(x,t,t/a)u- {a 2 (t/a)+a 3 (t/a)e' (3l ^ x ^}sat(x 2 ) V ' 

for any g satisfying our assumption Aq for some c a > and with a\(s) = s 2 /2, provided the 
constants c a > and a > are sufficiently large. Similar extensions can be made for the other 
examples we considered above. 



8 Conclusion 



We provided general conditions under which slowly time-varying systems are uniformly globally 
asymptotically stable and input-to-state stable with respect to general perturbations, thus extend- 
ing [PA2, SI] to situations where the corresponding frozen dynamics are not necessarily exponen- 
tially stable. Moreover, we provided new methods for constructing explicit closed form strict ISS 
Lyapunov functions for slowly time-varying control systems in terms of a suitable family of gen- 
eralized Lyapunov like functions for the frozen dynamics. This is significant because Lyapunov 
functions play essential roles in robustness analysis and controller design. 

We conjecture that our work can be extended to systems that are subjected to actuator or 
measurement errors, or which are components of larger controlled hybrid dynamical systems. It 
would also be of interest to extend our work to slowly time-varying systems with outputs and 
to construct corresponding input-to-output stable (IOS) Lyapunov functions; see |SWll ISW2] for 
further background on systems with outputs and [MMlj for some first results on constructing IOS 
Lyapunov functions for non-autonomous systems in terms of given nonstrict Lyapunov functions. 
We leave these extensions for future papers. 
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